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ON \ General discrete one-dimensional stochastic models to describe the transport of 

single molecules along coupled parallel lattices with period N are developed. Theo- 
retical analysis that allows to calculate explicitly the steady-state dynamic properties 
of single molecules, such as mean velocity V and dispersion D, is presented for N = 1 
and N = 2 models. For the systems with N > 2 exact analytic expressions for the 
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large-time dynamic properties are obtained in the limit of strong coupling between 
the lattices that leads to dynamic equilibrium between two parallel kinetic pathways. 
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I. INTRODUCTION 

Successful functioning of biological cells strongly depends on two classes of enzyme 
molecules, motor proteins and cytoskeleton proteins, that form the basis of the biologi- 
cal transport systems.— ^ Motor proteins such as kinesins, dyneins, myosins, RNA and DNA 
polymerases, etc., operate in cells by transforming the chemical energy of hydrolysis of ATP 
(or related compounds) into the mechanical work, although the mechanism of this process 
is still largely unknown.— The cytoskeleton proteins, such as microtubules and actin fila- 
ments, are rigid multifilament linear polymers that provide tracks for the motion of motor 
proteins . 1121 ? 14 The growth dynamics of these proteins strongly influences their biological 
functions.-i* 2 ^ One of the most striking properties of growing cytoskeleton proteins is the 
dynamic instability phenomena that observed both in vitro and in vivo in the microtubules 3 
and in ParM proteins,— closely related homologs of the actin filaments. 

Recent experimental advances have allowed to determine the dynamics of motor pro- 
teins and the growth of cytoskeleton proteins with a single-molecule precision at different 
conditions.-^^iLSiSiili These experiments suggest that the complex biochemical transitions 
and intermediate states in the motor proteins and cytoskeleton proteins influence their dy- 
namic properties and functions. In order to understand biophysical properties of these 
proteins all biochemical pathways and transitions should be taken into account. 

One of the most successful approaches to describe the dynamics of motor proteins 
and growth of cytoskeleton proteins is a method of multi-state discrete chemical kinetic 
or stochastic models . 1111211 ? 11411 ^^ 1 ^ 1 ^ 1 ^ 21 ? According to this approach, the motor protein 
molecule or the tip of the cytoskeleton filament moves between discrete biochemical states, 
and these transitions are governed by a set of stochastic rates. These biochemical pathways 
are periodic due to the periodicity in the structure of the cytoskeleton proteins. This method 
allows to obtain exact analytic expressions for the dynamic properties, such as mean velocity 
V and dispersion D, in terms of the transition rates for systems with arbitrary number of 
states in the periods . 12il3il4i] ^ Another advantage of the discrete stochastic models is the 
ability to describe the complexity of the underlying biochemical pathways, something that 
cannot be done by other theoretical methods . 12i2 ° Currently, discrete stochastic models are 
developed, analyzed and applied for the simple sequential single-chain pathways^^^, for 
the biochemical pathways with irreversible det achment o 13 1 16 , for the lattices with branched 



3 



states,— for the parallel-chain pathways , 13 i 15 and for the systems with general waiting-time 
distributions for transitions between the states.— 

In order to describe better the complex biological transport processes, in this work we 
extend the discrete chemical kinetic approach by considering a two-chain model with a direct 
coupling between parallel lattices at each site, as illustrated in Fig. 1. There are total 2N 
sites in each period of the system, equally divided between the chains. Each site corresponds 
to a specific biochemical state of the the motor protein or the cytoskeleton protein. The 
protein molecule can be found on the upper biochemical pathway (chain 0) or it can move 
along the lower pathway (chain 1): see Fig. 1. The particle at state j (j = 0,1, ■ ■ ■ , N — 1) 
on the chain can make a forward step with the rate Uj and a backward step with the rate 
Wj — see Fig. 1. It also can make a vertical transition to the state j on the chain 1 with 
the rate jj. Similarly, the particle at state j on the chain 1 can diffuse forward (backward) 
with the rate <x,- (/3j), while the transition to the upper channel is given by the rate 5j, as 
shown in Fig. 1. The distance between identical states in different periods is equal to d. 
This model can be used, for example, to analyze the dynamics of RNA polymerases^ or it 
can describe the dynamic instability of microtubules and ParM proteins.-*^ 

Our analysis of the periodic coupled parallel-chain discrete stochastic models is based on 
the Derrida's method— that allows to obtain explicit analytic expressions for the stationary- 
state drift velocity 

V = \ '({//,. «•,.'> ,. )' ; .- r ^ J ) = lim ^-{x(t)), (1) 

and dispersion (or diffusion constant) 

D = DUu^Wj^j^j^^Sj}) = ~ lim | [(x 2 (t)) - (x(t)) 2 } , (2) 

z t^oo at 

where x(t) represents the spatial coordinate of the motor protein or of the tip of the cy- 
toskeleton protein. 

The paper is organized as follows. The results for discrete parallel-chain models in the 
limit of strong coupling between the pathways are outlined in Section II, while general 
explicit formulas for dynamic properties for N — 1 and N = 2 models are discussed in 
Section III. Section IV summarizes and concludes are analysis. The details of mathematical 
calculations are given in Appendix. 
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II. TRANSPORT ALONG PARALLEL CHANNELS IN THE LIMIT OF 

STRONG COUPLING 

Consider the general periodic coupled two-chain discrete stochastic model as shown in 
Fig. 1. Let us define J(j, 0) and J(j, 1) as probability density currents between the states 
j and j + 1 in the channel and 1, correspondingly. Also J{j,v) is a vertical probability 
density current from the state j in the channel to the state j in the channel 1, assuming 
that this current is positive in the downward direction (see Fig. 1). The stationary-state 
conditions and the periodicity of the system require that J(j, 0) + J(j, 1) to be a constant 
for any j = 0,1, ■ ■ ■ , N — 1. Then it leads to the following condition on the vertical currents, 

N-l 

E j c^) = - ( 3 ) 

3=0 

Generally, each vertical current is not equal to zero, and to calculate the dynamic prop- 
erties of the system analytically is mathematically very hard, except for N = 1 and N = 2 
models as will be shown below. However, when the vertical transitions rates 7j and 5j are 
much larger than other transition rates, one expects that in the limit of large times there 
will be a "dynamic" equilibrium between two channels, i.e., J(j,v) = at each site. We call 
this a strong coupling limit. In this shown in Appendix, it is possible to obtain exact 

analytic expressions for the dynamic properties of the system. Specifically, the equation for 
velocity is written as 



(4) 




and the expression for dispersion is 



d 2 I I 

D = A. 



N 2 



N-l N-l N-l , \ 

N W, + A s 3 + 1) 1 + ^ r j+l+1 
1^ 1^ 7^ V tj+i+i/ 



3=0 3=0 



where A = NV/d, and the modified transition rates are 

7j" .7. , ^ 



= Uj + yaj, vjj = Wj + -ffij. (6) 

°3 °3 
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Also, the functions rj and Sj are defined as 



and 




Comparing these results with the velocity and dispersion for single-channel periodic dis- 
crete models , 12 ! 21 it is interesting to note that the dynamic properties of the two-channel 
coupled model can be viewed SIS db 9bll effective motion along the single pathway with the 
modified transition rates. This is due to the dynamic equilibrium between two parallel 
pathways in the strong coupling limit. 



III. DYNAMICS FOR N = 1 AND N = 2 MODELS WITH GENERAL 

COUPLING 

For the periodic two-chain coupled discrete stochastic models with N — 1 the dynamics 
can be calculated explicitly for all parameters. Because of Eq. Q the vertical current is 
always equal to zero at all sites, and the dynamics is the same as in the strong coupling 
limit. Then the expressions for the mean velocity and dispersion can be easily obtained from 
Eqs. © and ©: 

7 + 5 

D _ d 2 1 (u + w) + 5(a + (3) 

2 1 +5 ' 1 ' 

It can be seen that both mean velocity and dispersion consist of two terms, corresponding 

to the motion along the pathways and 1. Here and are the probabilities to find 

the particle on the channel and 1, respectively. 

Now let us consider N = 2 periodic two-chain discrete stochastic models. In this case 

the dynamic properties can be calculated by solving explicitly a set of Master equations (see 

Appendix). Then the exact expression for the mean velocity can be written as follows, 

v = Ti {(«o«i - wowi) [(«o + Po)5i + (on + Pi)S + SoSi] 
+ (a «i - AA)[(«o + w )7i + {u\ + wi)7o + 7o7i] 

+(ii ai - ioq/3i)£o7i + ("o^i - A^ihiA} , (11) 
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where 

Q = (u + w + 7o)[(«o + Po)8i + (ati + ^1)71] 
+(ui + wi + 7i)[(«o + Po)lo + («i + Pi)S ] 
+(«o + p + 6 )[(u + 100)71 + («i + w^Si] 

+(«! + A + (Ji)[(tto + w )5 + (mi + 101)70]. (12) 

The corresponding expression for dispersion is quite bulky and it is presented in the Ap- 
pendix. 

It is interesting to compare these formulas with the results from the strong coupling limit. 
The mean velocity is given by a simpler equation, 

(uq$o + Qo7o)(mi(u + ornj - (w 5 + /3o7o)(w £o + /3q7q) 

(uq + w )5 + (a + /3 )7o)(7i + <*i) + (ui + ^1)^1 + («i + A)7i)(7o + ^o) ' 

This formula directly corresponds to Eq. (j3J). 

In order to illustrate the effect of the coupling between the pathways we represent the 

vertical transition rates in the following form 

Ti(e.P) = 7°ep. ^-(e, p) = - p). (14) 

We take 7° and 5° as fixed parameters and vary only the parameters e (e > 0) and p 
(0 < p < 1). The coefficient e gives a measure of the coupling between the channels. When 
e > 1 we reach the strong coupling limit, while for small values of e the coupling is rather 
weak. The coefficient p reflects the direction of the dominating vertical current. For p ~ 
the current is mainly upward from the chain 1 into the chain 0, while for p ~ 1 the vertical 
current changes the direction (see Fig. 1). The mean velocities for N = 2 model with 
different sets of parameters are presented in Fig. 2. It can be seen that the mean velocities 
do not depend strongly on the coupling between the pathways. Thus the expression for the 
mean velocity generally can be well approximated by the strong coupling result. 

The situation is different for the diffusion constant. Dispersions for different parameters 
are plotted in Fig. 3. It is interesting to note that for a given set of the transition rates 
dispersion reaches a maximum for low values of e and intermediate values of the parameter 
p. When the coupling between the channels is weak (e = 1) the particle spends most of the 
time by diffusing in one of the channels before transferring to another one, and dispersion is 
large. However, for strong coupling (e — > 00) the transitions between the vertical states are 
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very frequent. The particle moves a short distance along the channel before transferring to 
the other chain, and, as a result, dispersion is low. The situation for low couplings between 
the channels is very similar to the picture of dynamic instability in microtubules.—^ Note also 
that this behavior cannot be observed in iV = 1 models. This suggests that intermediate 
biochemical states might play a critical role for understanding dynamic instability in the 
cytoskeleton proteins. 



IV. SUMMARY AND CONCLUSIONS 



Parallel-chain periodic discrete stochastic models with coupling between the channels are 
introduced and studied theoretically. The exact analytic expressions for the mean velocity 
and dispersion are found for the models with arbitrary period size in the limit of strong 
coupling. This corresponds to the dynamic equilibrium between the pathways and zero 
currents in the vertical directions at each state. For the two-chain discrete models with 
N — 1 and N = 2 periods the dynamic properties are calculated explicitly for general 
conditions. It is found that for weak coupling between the channels dispersion is maximal, 
while in the dynamic equilibrium regime it is significantly lower. At the same time, the effect 
of the coupling on the mean velocities is much smaller. Our analysis can be extended to more 
general parallel-chain discrete stochastic models with more than two different pathways. In 
addition, it will interesting to to study the dynamics of these models under the influence 
of external forces.— It is also suggested that these models can be used to investigate the 
transport of motor proteins and growth dynamics of cytoskeleton proteins. 
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Appendix 

Let us introduce the probabilities Po(l;t) and Pi(l;t) of finding the particle at time t at 
the position I = j + Nk (j = 0, ...N — 1, k are integers) on the lattice and 1, respectively. 
These probabilities satisfy the following Master equations, 

dP (j + Nk,t) 



dt 



Uj-iPoO" -l + Nk,t)+ w j+1 P (j + 1 + Nk,t) + SjP^j + Nk, t) 



-{ Uj + Wj + ij)P (j + Nk, t), (Al) 

dPi(j + Nk,t) 



dt 



aj-iPiO" — I + Nk, t) + p j+1 Px(j + 1 + Nk, t) + 7i P (j + Nk, t) 

-{aLj + Pj + 5j)Px(j + Wfc, t). (A2) 



The conservation of probability requires that 

+oo N-l +oo N-l 

k=—oo 3=0 l=—oo j=0 

at all times. 

Following the idea of Derrida,— we define two sets of auxiliary functions for each j = 
0,1,- ■■ ,N-1 audi = 0,1. 

+oo 

Bi(j,t)= Pi(j + Nk,t). (A4) 

k=— oo 

Ci(j, *) = X) V + ^^0' + *)■ (A5) 

fc=— oo 

Note that the conservation of probability yields 

N-l N-l 

J2B (j,t) + Y,B 1 (j,t) = l. (A6) 

3=0 3=0 

Then from the Master equations ()A1|) and ()A2|) we derive 
dB (j,t) 



dt 

dB^t) 
dt 



uj-xBqU -l,t) + w J+1 B (j + l,t) + SjBxij, t) - ( Uj + Wj + ^)B (j, t), (A7) 



a^xB^j -l,t)+ P j+1 Bx(j + l,t)+ ijBoij, t) - (a, + fa + S^B^j, t). (A8) 



Similar arguments can be used to describe the functions C (j,t) and C x (j,t): 
dC (j,t) 



dt 



Uj-iC (j - 1, t) + w j+1 C (j + l,t) + SjCxij, t) - (uj + wj + 7j)Co(j, t) 



+u j - 1 B {j - 1, t) - w j+1 B (j + l,t), (A9) 

= -l,t) + /3 i+1 C 1 (j + 1, t) + 7j C (j, t) - (a 3 + (3, + SJCM, t) 

+a j -. 1 B 1 {j - l,t) - Pj+ x B x {j + l,t). (A10) 

Again using Derrida's method, 21 we assume the following stationary-state behavior 
Bi(j, t) - bi(j), Ci(j, t) - Oi0)t + W) (» = 0, 1; j = 0, 1, • • • , N - 1). (All) 
At steady state we have dBi(j,t)/dt = 0, and Eqs. (|A7J) and (|A8|) yield 

= u^-i&oCj - 1) + Wj+ibo(j + 1) + (j) - + Wj + Jj)b (j), (A12) 

= otj-i&iO' - 1) + Pj+Mj + 1) + 7^o(j) - (a, + P 3 + SjMj). (A13) 

To determine the coefficients at(j) and Tj(j), Eq. (|A11J) is substituted into the asymptotic 
expressions (jA9|) and (jAlOj) . producing 

= Wy_ia (j - 1) + u;j+ia (j + 1) + ^ai(j) - (itj + Wj + TjHO'), (A14) 

= aj-iaiO' - 1) + Pj+iai(j + 1) + 7i«o(j) - Oj + #7 + ( A15 ) 
The coefficients Ti(j) satisfy the following equations 

«o(j) = Uj-iT Q (j - 1) + w j+1 T (j + 1) + 5,-TiO') - (uj + wj + ij)T (i) 

+Uj-ib (j - 1) - tu i+1 & (j + !): ( A16 ) 

= "i-i^iO' - 1) + Pj+iTiU + 1) + 7i^o(j) - («i + /3j + ^-)Ti(j) 

+a i _ 1 T 1 (j - 1) - ^+1^0' + 1). (A17) 

The strong coupling limit describes the situation when the rates Sj and 7j are much larger 
than other transition rates. In this case we have 

B^j^-^B^t), C x {j,t)~^C Q tj,t). (A18) 

°3 °j 
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Then, by introducing the modified transition rates 

Uj = Uj + -jr-otj, Wj = Wj + -jf/3j, (A19) 
from Eqs. (lA"T2l) . ljXT5|) . and ijXTSj) we obtain 

= Uj-Mj - 1) + %+i6o(j + 1) - (uj + %)6o(i)- (A20) 
At the same time, Eqs. fA"T4|) - (jXTTj) . and fA"T8|) lead to 

= Uj-idj-i + w^+iOj+i — («j + Wj)dj. (A21) 

(1 + ^-)a(j) = Mj-iTj-i + w j+1 T i+ i - + u^T.,- 

- 1) - + !)■ (A22) 

Here for simplicity we put ao(j) = a j and To(j) = Tj. Following Derrida's solution,— we 
show that bo(j) = Qo r j with 



, N-l k . 

1 i+En^i < A23 > 



M J V L._1 4 _1 



fc=l i=l 

The coefficient Go can be found from the normalization condition (|A6|) and from Eq. (|A18|) . 

o = N _, - 1 - • (A24) 

E (i + lU 

i=o v J/ 

Then, comparing Eq. (|A21J) and Eq. (|A20J) . we conclude that aj = Abi(j). The coefficient 
A can be calculated by summing up Eqs. (jA22|) for j = 0, 1, • • • , N — 1, 

N-l , v N-l 

A = E ( 1 + t) % = ~ ^ >o(i) ' (A25) 

that yields 

tiV-l 
i- n | 

- = ^F— r^- (A26) 



3=0 V J/ 

To solve Eqs. (jA22|) we define 

Vj = *i+i T i+i - % T 3- (A27) 
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Then Eqs. (jMj) 

can be rewritten as 
yj - = ^1 + j-J a j + Wj+ibod + 1) - Uj-ib (j - 1). 
The solution of this equation is 



(A28) 



Vj = U j 



& o(j) + ^E^ + 1 ) (i + ^ Uo(j- + i + i) + e 



(A29) 



where Bi is some unknown constant. From Eq. (jA27|) we obtain 

N-l k 



T; 



j+i 



k=l i=l 



Uj+i 



(A30) 



Now we can can calculate the mean velocity using the expression for the mean position 
of the particle 



+oo N-l N-l 

< x(t) >= | E E U + Nk)[PoU + Nk,t) + P 1 (j + Nk,t)) = | £ [C (j,t) + 

fc= — OO jr=0 j=0 

(A31) 



2V 



e 1 (i+|)c a,t). 



The velocity is given by the following equation 

N-l 



J=0 v J ' 7=0 v 



where the function A is given in Eq. (jA26j) . 

Similar calculations can be performed to determine dispersion. Starting from 

+oo N-l 

<x 2 (t)>= ^(j + iVA;) 2 [Po(j + iVfc,t)+Pi(j + iVA;,t)], 

k=— oo j=0 



U) = jr A > ( A32 ) 



(A33) 



and utilizing Eqs. flXTJ), (jA"2j) and ()A19j) . we obtain 

N-l 



d 



N-l 



lim — < x 2 (t) >= 2 Vi 



J'=0 



)C (j, t) + + Wj)B {j, t). 

3=0 



(A34) 



Then, the diffusion constant can be derived from the definition (J2j) 

"JV-l , N-l N-l 



d 2 



JV— 1 ^ i V — 1 iV — 1 y \ 

- ^ + 2 + - A E ( 1 + 1) T ^ 

i=o j=0 j=0 ^ 3/ 



(A35) 
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By substituting the expressions for 7} [using Eqs. (jA30|) and (|A29|) ] into Eq. ()A35|) . the 
final expression for dispersion, given in Eq. (jSJ), is obtained. The unknown constant 0i 
cancels out in the final equation for dispersion. 

For N = 2 models with arbitrary coupling between the channels the dynamic properties 
can be obtained by solving directly Eqs. (|A12j) and (|A13j) . The solutions are 

(ui + tui)(a + 0o + 5 )5i + (ai + (3i){ui + w x + 7i)£ 



6o(0) 
6i(0) = 



n 

(u + w )(ai + (3i + Si)S + (a + f3 )(u Q + w + j )5i 

n 

(ai + Pi)(uq + w + 70)71 + (m + Wi)(ai + Pi + S^p 

n 



6i(1) = . (A36) 

The expressions for di(J) follow from the relation di(J) = Abi(j). Similarly, the general 
solutions of Eqs. (jA16|) and (jA17|) yield the explicit expressions for Tj(j'). These expressions 
are not shown here because they are quite bulky. Finally, the formula for the velocity is 
given in Eq. (jllj) . while for dispersion we have 

d 2 

D = j {K - ^o)To(O) + ( Ul - Wi)To(l) + (a - A))2i(0) + (a a - ^^(l) 

-A(T (0) + T (1)+T 1 (0)+T 1 (1)) 
+ w )b (0) + (ui + iwi)6o(l) + (ao + Po)h(0) + (on + (3i)h(l)\ . (A37) 
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Figure Captions: 

Fig. 1 General kinetic scheme for the two-chain periodic stochastic model. Both channels 
have N discrete states per period. The particle on the upper chain moves forward (back- 
ward) with the rate Uj (wj), while on the lower chain 1 the corresponding rates are aij and 
(3j ( with j — 0, 1, • • • , N — 1). The vertical transitions between the pathways are given by 
7j and Sj. 

Fig. 2 Comparison of mean velocities for N = 2 coupled parallel-chain discrete stochastic 
model as a function of the parameter < p < 1 for different values of the parameter e (see 
text for explanations). Solid line corresponds to e — 1, dashed line is for e = 20, while 
the dotted line describes the strong coupling limit (e — > oo). The values of the horizontal 
transition rates used for calculations are u = (3 = 5, w = a = 1, u± — (3± — 10, 
w 1 — a± — 2. For transitions between the chains we used the following parameters: 7$ = 1, 
^0 = 7i = 2, and <5° = 0.6. For all calculations it was assumed that d—1. 

Fig. 3 Comparison of dispersions for N = 2 coupled parallel-chain discrete stochastic 
model as a function of the parameter < p < 1 for different values of the parameter e (see 
text for explanations). Solid line corresponds to e — 1, dashed line is for e = 20, while the 
dotted line describes the strong coupling limit (e — > 00). All parameters used for calculations 
are the same as in Fig. 2. 
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